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ABSTRACT 


Equations are developed for synthesizing the radar cross section (RCS) of a 
planar resistive sheet. The equations are verified by first calculating the RCS of 
the sheet for an assumed resistivity function. The computed RCS is used in the 
synthesis equations, and the resistivity distribution is shown to be in agreement 
with the original assumed distribution. 

Several synthesis parameters are examined and their effect on the accuracy 
and computing time are discussed. Results are shown for both monostatic and 


bistatic RCS with arbitrary incident wave polarization. 
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I. INTRODUCTION 


The design of radar-stealthy platforms has become an 
important engineering problem over the last several decades. 
Although the principles of radar stealth are well known, radar 
cross section (RCS) prediction and reduction are still 
difficult tasks. The development of increasingly sophisticated | 
detection systems threatens to reduce the mission 
effectiveness of many types of weapons platforms. Much effort 
is now spent on reducing radar detectability by reducing the 
RCS of the platform. 

' Four basic techniques are employed in radar cross section 


reduction (RCSR): 


1. shaping, 
2. surface material selection, 
3. active cancellation, and 


4. passive cancellation. 


Fach method has advantages and disadvantages and the final 
selection is a compromise. Shaping is one of the most 
important ways of controlling the radar echo, but usually RCSR 
at one viewing angle is accompanied by an enhancement at 
others. Radar absorbing materials (RAM) can also be used, and 
the RCS reduction for this case is achieved by the dissipation 
of energy inside of the material. The platform performance 
(e.g., aerodynamic efficiency) generally decreases due to 
added weight. Surface maintenance requirements are also more 
demanding than those for conventional materials. Active and 
passive cancellation are limited methods and can only be 
applied to narrow frequency bands and spatial regions. 

When choosing materials it is necessary to establish a 
relationship between the electrical characteristics of the 
material and the target’s scattered field. The electrical 


properties of a material are usually expressed in terms of the 
complex permittivity and permeability (e€ and p respectively). 
Finding the scattered field when €, pw and the shape of the 
target are known, is straightforward, although not always 
easy. RCS synthesis, the inverse problem (i.e., finding ¢ and 
wp for a given RCS pattern and body shape) is more difficult 
because, in general, a complicated set of coupled integral 
equations must be solved. Because of the complexity of the 
synthesis procedure, RCS reduction using RAM has been based on 
intuition and experience. 

An RCS synthesis method for an impedance surface can be 
based on the solution of the scattering integral using a 
method of moments (MM) technique. The (unknown) surface 
impedance is expanded in a series of basis functions with 
unknown coefficients. In the most general Synthesis case, both 
electric and magnetic currents must be considered. When the 
Surface 1S planar, an expression for the currents can be 
obtained by applying an approximate set of boundary conditions 
at the surface [Ref. 1]. In general, a coupled set of integral 
equations must be solved [Ref. 2]. The solution yields the 
surface impedance at every point on the surface of the body 
for a specified scattered field (or RCS) pattern. Once the 
surface impedance is known, the electrical properties of the 
material can be determined under certain conditions via the 
refractive index of the material [Ref. 3]. 

A resistive sheet is a special case of an impedance 
surface since it does not support magnetic currents. It is of 
great importance for RCS reduction purposes because it is used 
to control travelling wave effects. In this thesis, the 
Synthesis procedure described above is applied to a resistive 
sheet, in which case the integral equations Simplify 
considerably. 

Chapter IIT is a review of the concepts and formulas 


required in the formulation and solution of the problem. The 











reflection and transmission coefficients for a thin resistive 
surface are also derived from the boundary conditions. In 
Chapter III the appropriate equations are derived and 
expressed in matrix form for the monostatic and bistatic 
cases. Chapter IV deals with the computer implementation of 
these equations and the evaluation of the obtained data. 


Finally, Chapter V concludes with a discussion on the 


benefits, concerns, and recommendations for future work. 








II. THEORETICAL BACKGROUND 


In this chapter, a few concepts and equations used 


throughout the thesis are introduced. 
A. RESISTIVE, CONDUCTIVE AND IMPEDANCE SURFACES 


Thin layers of lossy materials are of special interest in 
RCSR design as they have lower RCS than the corresponding 
perfect electric conductors [Ref. 3]. A mathematical model of 
such a layer is an electrically resistive sheet. It can be 
visualized as a thin sheet of highly conducting material whose 
permeability is that of the surrounding medium (free space). 
Theoretically, it is composed of an electric surface current 
of infinitesimal thickness whose total strength is 
proportional to the tangential electric field at its surface. 
The sheet is characterized by a jump discontinuity in the 
tangential components of the magnetic field across the 
surface, but no discontinuity in the tangential electric 
field. Therefore, a resistive sheet does not support magnetic 
currents. Its properties are completely specified by its 
resistivity, R, in ohms per square. This unit is derived from 
the basic definition of the resistance of a rectangular cell 
with length 1, and cross section, S. For the special case of 


length equal to width (1l=w) 


l 1 a Q) 
R= = _ = — 1 
© 0S atw oat co! 
where, o is the conductivity of the material in 


(Siemens/meter) and t is the thickness of the cell. 

The boundary conditions at the two sides of the sheet 
follow from the usual boundary conditions of electromagnetics 
[Ref. 2]. For the limiting case of t approaching zero, 0 


increases in such a manner that R, is finite in the limit. 


However, when R,=0 the sheet is perfectly conducting, and 





when R,=° it is perfectly transparent. Thus, the boundary 


conditions for the resistive sheet are [Ref. 4] 


where the signs refer to the upper (positive) and lower 
(negative) faces of the sheet, and J, is the total surface 
electric current supported by the surface. 

A general vector form of the boundary conditions are 


[Ref. 1] 


AxE(+) -AxE(-) =0 (2) 

AxH(+) -AxH(-) =F, (3) 
with 

Axt{AxE(+)}=-R,J, (4) 


where f is the unit vector normal to the sheet and directed 


towards the positive (plus) side of the sheet and J, is the 


total electric current supported. 

The electromagnetic dual of the resistive sheet is the 
magnetically conductive sheet having a conductivity R, in 
S/sq. Such a sheet has a jump discontinuity in the tangential 
electric field but none in the magnetic field, and Supports 
only magnetic currents. It simulates a thin layer whose 
permeability differs from that of the surrounding medium (it 


would be difficult to realize this type of sheet). When R,,=0 


the sheet acts like a "perfect ferrite" with infinite 


permeability and when R,=~ it no longer exists. The general 


boundary conditions for such a sheet are 


with 
Ax{AxH(+)}= -R Jy, 


where J, is the magnetic current supported and, the convention 


followed for the signs is the same as in the resistive sheet 
case. 

Note that an impedance surface is a simulation of a thin 
layer whose permittivity and permeability both differ from the 
surrounding medium. Therefore, a combined sheet consisting of 
resistive and conductive ones is equivalent to an impedance 
surface. Although the two sheets are, in general, coupled, 
with each affecting the scattering from the other, decoupling 
occurs when the sheets lie in a plane. It is therefore 
sufficient to solve the problem for the simpler case of the 
resistive sheet and use the duality principle to get the 
solution for the other. The solution of the combination sheet 


can follow by adding the two solutions [Ref. 1]. 
B. REFLECTION AND TRANSMISSION 


The reflection and transmission coefficients for a 
resistive film can be derived from the boundary conditions in 
Equations (2), (3) and (4) and the geometry shown in Figure 1. 
The subscripts i, s and t refer to incident, scattered (or 
reflected) and transmitted, respectively. The two principal 
polarizations (parallel and perpendicular) are investigated 
separately. Specular reflection, and transmission at the angle 


of incidence is assumed, because the sheet is infinitely thin. 





> — 
Bi= By 
Figure 1: Reflection and transmission from 
a resistive sheet. 


1. Parallel Polarization 
For parallel polarization, the electric field of the 


incident wave lies in the plane of incidence 
Bi=O6 Rie 7B: | (5) 


Without loss of generality, it can be assumed that the plane 


of incidence is the x-z plane. Figure 2 shows a cut through 


the @=0 plane. 





Figure 2: A cut of Figure 1 in the plane of 
incidence. 








The incident fields are 


Bizk [-Rcos(0,) +2sin(8,) ] e 78 msn Os) 2008 Wa)! (6) 


and 
=i »~ Ho _-Jp [-xsin(6,) -zcos (6,) ] 
H*=p—e 2 : (7) 
Oo 
where p= 2 is the wave number. Denoting fT, as the 


transmission coefficient for parallel polarization, the 


transmitted fields are 


fitz PE. [-2cos (8,) +8sin (6,) ] e JP [-xsin(@,) -zcos (6,) ] (8) 
and 
= ¢ Ey -7p [-xsin(6,) -zcos (6;)] 
H ae eC = + : (9) 
oO 
Finally the scattered fields are 
E°=TLE, [-%cos (6,) -2sin(6,;) ] eo dBi eet Os) S208 0a). (10) 
and 
pe va -jB [-xsin (6;) +zcos (6;) J 
H°=-yV,—e : (11) 


No 


where JT, is the reflection coefficient for parallel 


polarization. Using the notation introduced in Section A and 


applying the boundary conditions on the electric field 


intensity at Z=0 





(12) 


Similarly for the boundary conditions on the magnetic field 


intensity 
H(+) = (H?+H*)|,, 
(13) 
BS )e H* ines 
Inserting Equations (6), (8), (10), and (12) into (2) gives 
Lo Lal: (14) 
Using Equations (7), (9), (11) and (13) into (3) gives 
E . . 
= _T _ -jB [-xsin(6,)] 
J = aut Le oe : (15) 
Using Equation (4) and Equation (15) 
R 
=— = (1-I,-T,) . (16) 


T,,cos (6,) ; 


Solving Equations (14) and (16) for the reflection coefficient 


yields 
_ -79,cos (6,) (17) 
P  2R,+n,cos (6;) 
and for the transmission coefficient 
(18) 


a i 
P 2R,+n,cos (6;) 
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2. Perpendicular (Normal) Polarization 
For perpendicular polarization a similar procedure is 
followed for deriving the reflection and transmission 


coefficients. In general the incident field is of the form 


Bi=G ze TBH: | (19) 


Again, without loss of generality, the x-z plane is considered 
as the plane of incidence. Expressions for the incident waves 


are 


Bin OR e 7h here econ tal! (20) 


a ee os ~48 (-xe4 = 
Hi=- 2 [-2cos (0,) +28in(0,)] e JRE RRP nas (21) 


Oo 


Denoting T, and I, the transmission and reflection coefficients 
respectively for normal polarization, the scattered and 


transmitted fields are 


Be=-pr_ Bie IB mein (Os) ~zcos (8;) ] (22) 


= E ji - 7 -xS] ja ; 
H*=T,—* [-%cos (0,) +2sin(0,)] e FRE xsin Os) -ze0s 0) | (23) 


oO 


ee ie Da mana (24) 


and 


11 


H*=T,—° [cos (6,) +2sin(6,)] e978 ein.) +ze08 (8))1 | ggg) 


Nn 
oO 


E 
uf 


Applying Equations (20), (22), (24) and (12) to (2) gives 
Dad os (26) 
To obtain the surface current, Equations (21), (23), (25) and 


(13) are inserted into (2) yielding 


- E eee 
J_-9 —2 cos (8,) (Parr aa ye 2h seen! (27) 


oO 


Using Equation (27) in the boundary condition (Equation (3)) 


gives 
R 
T,-— cos (6,) (1-7,-Ir.) . (28) 


Solving Equations (26) and (28) for the reflection and 


transmission coefficients gives 


-n, 


"” 2R,cos (0;) +n, a 
and 
2R.cos (6,) (30) 


L = —  .. 
" 2kR,cos (6,) +n, 


C. RADIATION AND SCATTERING EQUATIONS 


In general, electric and magnetic currents are the 
sources of radiated or scattered fields. For the case of a 
resistive sheet the magnetic currents are zero. The formulas 
Simplify further when only the far zone fields are of 


interest. 


12 








The formulas for the far - field calculation are derived 
by Balanis [Ref. 5] and are given in this section for the 
reader’s convenience. They are specified for the case of 
currents in the x-y plane of a Cartesian coordinate system as 
shown in Figure 3. 


The spherical components of the electric field intensity 


are 
E.=0, 
1 e-iBbr 
B,- - LEE” NoNe» (31) 
ne 7p e-ibr 
Ey = ATT og» (32) 
where 
Ng=| f [J,cos (8) cos () +J,cos (86) sin (9) ] et Bricos¥ gy/ dy! (33) 
& 
Ny= ff _{-J,sin() +J,cos () ] eIBr’cos¥ gy! dy! (34) 
and 


r’‘cosw =x’sin (8) cos (@) +ty’sin (6) sin(@) =x/uty’veg. (35) 


For the observation point appearing in Equation (35) the 
direction cosines are defined by 


u=sin(6) cos(@) , 
v=sin(8) sin(®) , (36) 


w=cos (@) . 
In Chapter III the currents derived in Section B will be 


13 


used in the radiation integrals (Equations (31) and (32)) to 
obtain the scattered fields. The currents are functions of the 
surface resistivity through the reflection (or transmission) 
coefficients. Thus, for the synthesis problem E, and E, in 
Equations (31) and (32) are known but, the currents J, and J, 


(functions of R,) are unknown. 





<< |x ——___» 


Figure 3: Coordinate system for scattering 
from a planar rectangular plate. 


All the formulas assume that the current is not affected by 
the presence of the edge. This turns out to be a valid 
assumption for resistive surfaces of finite size because the 


edges are "softer" than perfectly conducting edges. 
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III. FORMULATION AND SOLUTION 


The scattering from a planar rectangular resistive sheet 
with surface resistivity R, is investigated for the monostatic 
and bistatic cases. Using Equation (4), an estimate for the 
current distribution is obtained, and the total scattered 
field is calculated. To verify the synthesis equations, the 
calculated field is used to initiate the synthesis procedure 
for the surface resistance. 

The geometry under consideration is still that shown in 
Figure 1. The sheet lies in the x-y plane and is illuminated 
by a plane wave whose electric field contains both parallel 


and perpendicular components 


Et = (6 Ey +6 Ey ) e@ iBit (37) 


where 
-j 8, -7=7B [x sin(6,) cos(o,;) ty sin(8,) sin(o,) +z cos (6;) ] 


=jP [x u;t+y V;+Z w,] . 


On the surface of the sheet (Zz=0 and F=f’) 
- iB, P=-7 8, -P’=7 B (x/u,ty’v,) =j7 Ph. (38) 


The formulation as well as the analysis of the problem is 
Gifferent for the bistatic and monostatic cases and, 


therefore, each is considered separately. 
A. BISTATIC CASE 
Equation (4) with f=2 becomes 
8x (Bx E(-)) =-RJ (39) 


Ss 


where 


15 





E(-) =E°= (OT, By +$T, Ey) eF8". (40) 
The cross products are 

2x6 =-2Rcos (8) sin() +¥cos (8) cos () 

2x =-£cos (b) -Psin(d) 


and the triple products 


Bx (2x6) =~ (cos (6) cos (db) +¥cos (6) sin () ) 
and 
2x (2x) =ksin(o) -¥cos(o) . 
Then 
2x (2x E(-)) ={-2[cos (6,) cos (o,) T,£y -~sin(o,) T, Ey] 


-¥ [cos (8;) sin(o;) T, Ey +cos (,) T, By ]}e7B*. 


Substituting in the explicit expressions for the transmission 
coefficients (Equations (18) and (30)), into (40), Equation 


(39) gives the components of the current on the surface of the 


sheet 
— cos (8,) cos(,) Ey sin(o,) cos (6,) E oi bh (41) 
‘a 2R,+n,cos (6,) 2R,cos (0,;) +1,) 
7 =9] £28 (6,) sin(o,) Ey , COS (p,) cos (6,) Es ei Bh. (42) 
y 2R,+7,cos (8,) 2R,cos (0;) +n, 


Now that the currents are known, the scattered fields can be 
calculated. Because of the complexity of the formulas, 
parallel and perpendicular polarizations are treated 
Separately. Besides, any arbitrary polarization can be 


decomposed into a sum of these two orthogonal components. 


16 





1. Parallel Polarization 
When the field is parallel polarized E,’ =0 and Equations 
(41) and (42) simplify to 


1 


= toi ph en eee Se 
J,,=2 Hy eF"cos (0,;) cos (;) 2R,*1,c08 (6,) ' (43) 
= 2 Fp eI BH eg Sia ae 
J,=2 Ee e7"cos (8,;) sin (g;) 2R.+n cos (6,) (44) 
Using equations (43) and (44) in (33) gives 
Ng = 2 Ea cos (8) cos (8,;) cos (o-;) I, (45) 


where a new quantity has been defined for convenience 


e Jk (h+g) J / 
=f, 2R,+n,cos (8, Zea cos Oi pee 


The resistivity is a function of position in general and, 
therefore, it cannot be taken out of the integral. In order to 
reduce the scattering equations to closed form expressions, it 
is assumed that the resistivity profile can be accurately 
represented by a step approximation. The resistivity function 
is represented by a series of basis functions (pulses) with 


unknown coefficients 
M N 
RAVES ye Set ey (47) 


where M and N are the number of patches in the y and x 
directions respectively. The unknown expansion coefficients 
are {r,,} and the expansion functions are two-dimensional steps 


or "pedestals" 
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13. xyes 
f. (x’, y’) = y me (48) 


0, else. 


The geometry is defined in Figure 4. The center of the mn® 


patch is at (X,,, Yan) and each patch has dimensions Ax by Ay. 


—» «—Ax 





Figure 4: Division of the sheet into 
rectangular subdomains. 


Applying the series representation (47) to (46) gives 


e JB (Arg) foe 
oe Tages 


7B (h+g) 
=f a ey 
2) > Lis LmntNecos (0;) 


e JB (h+g) 


ax’ dy". 


For the mi” term of the sum the integral is 


eJb (h+g) / J 
Soe eee ele dy/= 
If 2 ImntN,cos (6,) ened 


(49) 


1 3k (h+g) 
2 Ian tN cos (8 ay I. . eae 
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Evaluating the integral in (49) yields a closed form 


expression in terms of sinc functions. Thus J, becomes 


M N 
I,=AxAyS,S,))> a,nerr™ (50) 
m=1 n=1 
where 
_ sin[B (utu,;) Ax/2)] __, 6B (utu,) Ax 
ae CTT 7 ec a 2 
_ sin[B (vtv,) Ay/2) __, B (vtv,) Ay 
i eae f (52) 
_ 1 
Cnn 27.4) ,C0S\0y) 2H 
and 
P= Xmn (UtU;) t¥mn (VIV:) » (54) 


Using Equation (50) allows Equation (45) to be written as 


Ny =2Ee cos (8;) cos (@) cos (-o,) AxAyS,S, 


m=1 n=1 


Finally, the complete expression of the @-component of 
the scattered electric field for a parallel polarized incident 
field is 


Egi, = C, Ee cos (8,) cos (8) cos (@-9,) I, (55) 


where 


1 


iyecee ne, (56) 
The first subscript on E,’ denotes the polarization of the 
scattered field, whereas the second subscript denotes the 
polarization of the incident wave. 
By the definition of the radar cross section for a three- 


dimensional target 


G72 


aap 
o=lim[42r? ap (57) 
E? 


one obtains for this case 


_ (1B)? 
T 


Ogo [cos (6,) cos (6) cos (o-,) I,]? (58) 


To calculate the ¢-component of the scattered electric 
field, E,, Equation (32) is used along with Equations (34), 
(43) and (44). Following the same procedure 


Exo = -C, Es cos (6,) sin(-,) I, (59) 


and the corresponding radar cross section is 


2 
O49 = Pe [cos (0,) sin(@-@,) I,]*. (60) 


2. Perpendicular Polarization 
When the incident field is normally polarized (EK, =0), 


the current components in Equations (41), (42) reduce to 


1 


J,=-2 EJ eib4sin(o.) cos (0. ee 
‘ 103) 19) 2R,cos (6;) +H, 


(61) 
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: (62) 


= i og iph 2 
J, =2 Eges@cos (b;) cos (0;) 2R,cos (0,) 1, 
Using Equations (61) and (62) in (33) gives 
Ng =2 Ey cos (6) cos (8,) sin(o-@,) TZ, (63) 


where 


}k (h+g) 
r= ff. sinless a ay 
2R,cos (8,) +n, 


Following the procedure used for parallel polarization one 


obtains 
M oN . 
T,=AxAyS,S, ¥ be (64) 
m=1 n=1 
with 
b ee eee (65) 


ue 27 .COSs.0,) 4 


Finally, using Equation (63) in (31) gives the 6-polarized 


component of the scattered field 
Egg = Cy Eg COs (6) cos (6,) sin(o-@,) I,. (66) 


The corresponding radar cross section is then 
_ (noB)? 9 = 2 (67) 
S99 = —_— [cos (8) cos (0,;) sin(o-@,) Z,]°. 


To determine the ¢?-component of the scattered field, Equation 


(34) is used, yielding 
Ny=2 Ey cos (6,) cos (o-@;) I, 
and then 
Eg, = C,E§ cos (8;) cos (-,) I, . (68) 


The corresponding radar cross section is 
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[cos (8,) cos (o-@,) I,]*. (69) 


(1B)? 
a 


3. Arbitrary Polarization 

In general, the incident wave is composed of both 
parallel and perpendicular components. Thus it is necessary to 
combine the results of Subsections 1 and 2. The total @ and ¢ 


components of the scattered field intensity are given by 


Ey’ = Ego + Eg 
| | (70) 
=C,cos (6) cos (6;) [Eg cos (o-,) I,+£ysin(o-,;) I] 


and 
Ss _ Ss S 
Ey = Ego t Ege 


(71) 
=-C,cos (0,;) [Ey sin (o-@,) I,-Eg cos (O-O,) I,] . 


At this point it is possible to cast Equations (70) and 
(71) in matrix form 














& 
Ee Z Age A9¢ S; (72) 
where, uSing Equations (50) and (64) 
Aga = Ep C,AxAycos (6,) cos (6) cos (o-6,;) $5, (73) 
Agg = Eg C,AxAycos (6,;) cos (8) sin($-9,) 5,5, (74) 
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Ago =~ Ep C,AxAycos (8;) sin (o-$,) SyS, (75) 


and 


Agg = EG c,AxAycos (8,) cos ($-9;) 5,S,. (76) 


The terms S, and S, in the right hand side of Equation (72) 
represent the double summation terms in (50) and (64) for each 


scattered field point, namely 


M oN 
54 =). »> annem (77) 
m=1 n=1 
and 
M N 
Se). ye ee (78) 
m=1 n=1 


The resistivity is embedded in the terms S, and 5). 


B. MONOSTATIC CASE 


For the monostatic case (¢=¢, and 6=6@;) and Equations (70) 


and (71) reduce to 

Eg? = CEs cos? (6) I, (79) 
and 

Eg = C, Eg cos (8) I, } (80) 
respectively. Equation (53) can be written now as 


1 


a pl read 
mi 22r,n+N,cos (9) 


(81) 


and Equation (65) as 
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be Geran (82) 
Equations (79) and (80) cannot form a linear system of 
equations of the form of Equation (72) because the 
coefficients in vectors a and b (Equations (81) and (82)) 
change with angle. 
The crudest approximation that can be made to eliminate 
the @ dependence is to consider cos(@) +1. Equations (79) and 


(80) reduce to 

Eg =C, Ep I, (83) 
and 

Es = C, Eg I, (84) 
and Equations (81) and (82) to 


4 
a=) =. (85) 
2 LintNo 


The monostatic RCS is for both polarizations 


2 
a= f0B) r. (86) 
T 


Now it is possible to set up a linear system of 
equations. In fact, as a @-polarized incident electric field 
creates only @-polarized scattered electric field, and 
Similarly for a ¢-polarized incident field, an equation of the 
form of (72) contains only one component of the scattered 
field. These equations were implemented in MATLAB and the code 
listings are given in the Appendix. 
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C. SYNTHESIS FORM OF THE SCATTERING EQUATIONS 


For RCS synthesis, the Scattered field components derived 
in the two previous sections are known quantities and the 
expansion coefficients {a,,} and {b,,}, must be determined. 
Expressing the scattered field equations in appropriate matrix 
form, it is possible to solve for the expansion coefficients 
using standard matrix techniques. To achieve this it is 
convenient to use a single index to count the patches rather 


than the double summation. Let 


MxN 


M oN 
> yen » c,e7P (87) 
kot 


where c., can be either a,, or b,,. 

There are MxN unknowns and therefore MxN equations are 
needed to solve for {a,,} OF {Dyn}. Because of the periodicity 
of the sinc functions (Equations (51) and (52)) in the 
direction cosines u and v, it is advantageous to work in 
direction cosine space. As shown in Figure 5, the region 
defined by -1ls u,v <1 is divided into the same number of 
subsections as there are patches on the sheet. The centers of 
these subsections (u,,v,) form a grid of points at which the 
scattered field components are specified, yielding MxN values 
for each field component. The area in the unit circle is the 
visible region while exterior points have no physical 


Significance. 
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Figure 5: Division of direction cosine 
space into subsections. 


The values of E, and &, are used to form a column vector 
of length MxN. Furthermore, at each of the patches, cos(@,), 
sin(¢,-¢,), cos(¢,-¢) can be expressed in terms of (u,,V,) 


according to the following relationships 
uptvyz=sin? (6,) , 


Vv 
a ees 


Inverting these formulas gives 


cos (8,) =/1- (ug+ve) , 


V 
gd, =arctan(—*). 
Ux 


(88) 


In the direction cosine space the terms cos(@), sin(¢-¢), 
cos ($-¢,) as well as the product S,S, can be written in the form 
a column vector of length MxN. Denoting column vectors with an 


underline 
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[x,] 7 COs (O,) , 
[*],  COF (od, d;) / 
(89) 
[4] =sin(,-6;) f 
[X,] =9,, (U)So2(Ve) 
We can further define the column vectors X, and X as 
MxN ; 
KJ, 7 wee (90) 
vel 
and 
MxN ‘BP. 
= Jd vy (Up, V, 
[S], 7 > De BED (91) 
vel 
The vectors X, and X, can also be written as 
X= Ba (92) 
and 
X= Bb (93) 


where a and b are column vectors of length MxN and follow 
from equations (53) and (65) respectively. The matrix B has 
dimensions (MxN) by (MxN) and is defined as 


LB] gy = 7 Bb (uxtus) +¥4 (verve) (94) 


where k,v=1,...,MxN. 
With the above matrices defined, Equation (70) can be 


written as 
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=C,cos (6,) X *X (Ey X *X + Ey X,*X,] (95) 


where * denotes an elementwise product and C,=C,AxAy. 


Equation (71) becomes 


dc i = aad 
Ey=-C,cos (8;) X, +(Ey X,*X,- Hy X, *X,). (96) 


s 
d 


For the monostatic case X,=X,=X and Equation (83) becomes 


==) 
Eo=C, Ey X,*X (97) 
and Equation (84) 
ES=C, Ey X,*X. (98) 


For this case u=u and v,=v, and the matrix B is defined now as 


[B] ey = O72 PL tate Me) (99) 


In all of these equations the goal is to solve for the vectors 
xX, & or X, which can easily be done using elementwise 


division. Once these vectors are known, it is possible to 


solve for a or b of Equations (92) and (93). 
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IV. COMPUTER IMPLEMENTATION AND DATA ANALYSIS 
A. GENERAL APPROACH 


In this chapter the surface resistivity synthesis 
procedure is simulated and computational accuracy issues are 
discussed. The simulation consists of computing the RCS (or 
scattered fields) of a rectangular sheet with known 
resistivity. Then the computed RCS is used in the synthesis 
equations to obtain the resistivity function. In principal, 
the synthesized distribution should be exactly the same as the 
original distribution used to compute the RCS. However, in 
practice, the agreement between the two depends on several 
Factors, such as the sampling method used to fill the vectors 
BE, and £,’ in Equations (95) through (98). 

From the field equations in matrix form, the scattered 
field components can be computed in any direction (@, ¢) when 
the resistivity function R, is defined. To form the vectors E; 
and E,*, MxN directions must be chosen to compute specific 
values. Once the vectors EH, and E,are defined, it is possible 
to solve these equations for the vectors X,, X,, and X. The 
final step in the synthesis procedure is the reconstruction of 
the expansion coefficients {r,,} from these vectors. 

One approach to choosing the observation directions is to 
define MxN points in direction cosine space as described in 
Chapter III-C and Figure 5. The square region defined by 
-l<u,vsl1 is divided into patches and their centers represent 
the chosen directions. Equally spaced points in u and v have 
been used in most calculations, but the spacings in these two 
principal axes can be different in general. The advantage of 
this approach is that equal sample spacings work well with 
scattering patterns that are periodic in u and v. 

At this point, two variations have been investigated. 
They are referred to as method (I) and method (II). In method 
(I) the observation points are restricted to be inside of the 


29 


unit circle (visible region), which makes physical sense. A 
disadvantage of this method is the necessity of forcing the 
number of patches on the sheet to be equal to the number of 
sample points in the visible part of direction cosine space. 
This condition becomes very restrictive as the spacing in u 
and v gets close. Consequently, if MxN points are to be 
contained in the unit circle, the -lsu,vsi region has to be 
divided into more than MxN patches, say N, by N,. However, 
despite this restriction, the method was successful, even for 
large sizes of the sheet. 

A modified version of this approach is method (II) in 
which the MxN direction cosine samples are spread over the 
complete range -lsu,vsl. For this case N,=N and N\=M. Note 
that the points outside the unit circle are not in the visible 
region and therefore one might expect that these values would 
not affect the synthesis result. However, an attempt to set 
the exterior points equal to zero results in a badly scaled 
matrix for the fields, and the matrix is almost always 
impossible to invert. The reason is that setting the exterior 
field values to zero is an attempt to synthesize a resistivity 
distribution for a discontinuous’ pattern. Realizable 
resistivity functions do not yield scattering patterns that 
abruptly go to zero. To avoid this problem the calculation can 
be extended to the exterior points using complex angles. In 


this case the expression 


cos (0,) =/1-(u,?+Vv;,*) (100) 


can take on complex values. At the corners |u,=|vJ=1 and, 


Jcos (0,))=1 


(101) 
arg(cos (@,)) = 


rola 


While the first method makes more physical sense, the 
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second is more flexible in choosing the number of far-field 
points (or patches), and is shown to work well over a large 
range of values of M and N. 

To verify the synthesis procedure, the scattered field is 
calculated using two different resistivity profiles: (1) a 
constant resistivity (R=3772) and (ii) a two-dimensional 


linear taper given by the equation 


R,(X,yY) =150([x|+]y]) Q. (102) 


For the field calculation, a large enough number of patches is 
used to sufficiently approximate the linear shape of Equation 
(102). The synthesis was tested using a few specific field 
points in the range of 9 to 256. This corresponds to M and N 
in the range of 3 to 16. Beyond that, the method still works 
well, but computational time increases significantly. Roundoff 
errors and numerical instability also affect the solution. 
Matrix B in Equations (92) or (93) has to be inverted and is, 
in general, an ill-conditioned matrix. It becomes singular as 
the number of patches increases. In such cases the use of the 
Moore-Penrose pseudoinverse matrix [Ref. 7] rather than the 
inverse, was found to give satisfactory results for larger 


numbers of patches. 


B. BISTATIC CASE 


Methods (I) and (II) are applied for the bistatic case, 
first for a pure TM, incident wave, and then for arbitrary 
incident wave polarization, which is the most general 
situation. For each resistivity function, the scattered fields 
are computed for sheets of sizes 1Ax1A and 5Ax5X. Data are 
presented for various numbers of patches and the two 
resistivity profiles. The plot of the expansion coefficients 
{r.,} (i.e., the resistivity profile used to calculate the 


scattered fields) is compared with the reconstructed 
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resistivity function. For small number of patches, the 
reconstruction is exact (one can not distinguish them when the 
plots are overlaid). But as the number of patches increases 
some numerical and sampling errors can be observed. In general 
both methods (I) and (II) work well in the entire range of 


Synthesis parameters investigated. 


1. Parallel Polarization 
Assume the incident wave is T, polarized (Ey =0) and that 


it arrives incident from the direction 6,=30° and $,=0°. As 


an example, the RCS pattern for a one wavelength square sheet 
1s shown in Figure 6. The field has been computed using method 
(IIT) and N,=N\=64. Various synthesis results are presented in 
Figures 7 through 17. Note that 


1. M and N are the number of subdivisions on the 
resistive sheet in the directions y and x (see Figure 4). 


2. N, and N, are the number of subdivisions of the region 
-lsu,Vvsl (see Figure 5). 


3. Lx and Ly are the dimensions of the sheet in x and y 
directions expressed in wavelengths. 


It should be noted that in method (I) N,xN,eMxN, while in 
method (II) N,=N and N,=M. In most cases the two methods gave 
identical results, thus the data will only be presented once. 

In Figures 7 and 8, only 16 field points are used with 
constant and linear taper resistivities. The initial and 
Synthesized distributions are identical. Figure 9 and 10 show 
that the synthesized linear taper resistivity has a maximum 
error of 1% when 256 points are used. The relative error for 
the constant resistivity case with the same number of points 
becomes 2% as shown in Figure 11. 

The RCS for the 5Ax5A sheet is shown in Figure 12. Again 
the two resistivity functions are considered. In Figures 13 
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(constant resistivity) and 14 (linear taper), only 16 field 
points are needed to recover the exact resistivity function. 
In Figures 15 through 17 it is seen that the error in the 
synthesized resistivity is also negligible when 256 points are 
used. It was found that for a sheet with dimensions kAxkaA, the 


choice of N,=Nj=k gives converged results. 
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Figure 6: Bistatic RCS of a resistive sheet with 
Lx=Ly=lambda, N,=N,=64, 6=30 deg, ¢=0 deg, R, is given by 
Equation (102). 
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Figure 7: Original and synthesized resistivity with method 
(I) are identical for this case. Lx=Ly=lambda, M=N=4, N,=4, 


Nv=5. 
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Figure 8: Reconstructed and original resistivity functions 
are identical. For this case Lx=Ly=lambda, M=N=4, N,=N,=4. 
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Synthesized and ‘original resistivity for 


Lx=Ly=lambda, N=M=16, N,=N\=16. The maximum error is 1%. 


Figure 9: 
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Figure 10: Difference of the synthesized and original 
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Figure 11: Synthesized and original resistivities. For 
this case, Lx=Ly=lambda, M=N=16 and N,=N,=16. The error is 


2%. 
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Figure 12: Scattered RCS from a resistive sheet. R, is 
given by Equation (102). For this case Lx= Ly=5*lambda, 
M=N=64, N,=N,=64, 6@;=30 deg, ¢,=0 deg. 
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Figure 13: Synthesized and original resistivity are 
jdentical. For this case, Lx=Ly=5*lambda, M=N=4, N,=4, 


N=5. 
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Figure 14: Synthesized and original resistivity are 
identical. For this case, Lx=Ly=5*lambda, M=N=4, N,=N,=4. 
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Figure 15: Synthesized and original resistivity. For this 
case, Lx=Ly=5*lambda, M=N=16, N,=N,=16. 
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2. Arbitrary Polarization 

For the general case where the incident field has both 
parallel and perpendicular components, both the E, and the E,° 
vectors given by Equations (95) and (96) depend on X, and X,. 
The vectors X, and X, are not independent and both implicitly 
contain the expansion coefficients {r,}. An expression of 
these equations in such a form that permits the solution for 
{r,,} Once the field vectors are known is desired. One approach 
is to specify the magnitude of the total scattered field as 


given by 


EB’ y [Fel “Ze 


but this leads to very complicated functions of a and b. 


Alternately, Equations (95) and (96) can be expressed as 


(103) 


Solving this system for X 


—_—a 


and X, yields 


_ [4 0g* £5 4g9*Ze | 


EAverers. : (104) 
i [404 A 


0° 494 *Ao0 | , 


and 


: A go* Hg A gg *Zy] 


[* >], [Ao *4 0 Age *4 60], 


(105) 
> 


Once X, or X is known, Equation (92) or (93) can be used to 
solve for the vectors a or b. The expansion coefficients can 
easily be determined using Equation (53) or (65). 

This procedure was applied to the resistivity functions 
of the previous section for sheets of dimensions 1x1 and 


5Ax5X. Method (II) was used for convenience and the results 
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are, for all practical purposes identical with the original 
resistivity. Figures 18 and 19 compare the original and 
synthesized resistivity profiles (linear taper) for the 1Ax1A 
Sheet and Figures 20 and 21 for the 5Ax5X and constant 
resistivity. In both cases 144 points have been used. For all 


the cases it was assumed that E, = E,' =1. 
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C. MONOSTATIC CASE 


The synthesis procedure can be applied to the monostatic 
case using Equation (97) or (98) for the scattered field. 
Alternatively, Equation (86) of the RCS can be used, which is 
the same for both polarizations. 

Methods (I) and (II) have been investigated for two 
resistive sheets of sizes 1Ax1A and 3Ax3A, and the two 
resistivity profiles considered previously; namely a constant 
of 3770, and a linear taper given by Equation (102). The two 
methods give identical results. 

Figure 22 compares the corresponding data for the 1Ax1X 
sheet with constant resistivity. In Figure 23 the difference 
between the synthesized and original distributions is plotted 


for a linear taper profile on a 4\x4X sheet. 
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Figure 22: Original and synthesized resistivities for the 
monostatic case. On this case, Lx=Ly=lambda, M=N=12, 


N,=N,=12. 
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Figure 23: Difference in the synthesized and original 
linear taper resistivity for the monostatic case. For this 
case, Lx=Ly=3*lambda, M=N=12, N,=N,=12. 








D. RELATION TO THE SAMPLING THEOREM 


In the previous two sections it was demonstrated that the 
synthesis procedure works well for both methods (I) and (ITI) 
over a wide range of synthesis parameters. These calculations 
have verified the synthesis equations and the computational 
approach to solving them. So far the question of how many 
field points are needed to uniquely determine an unknown 
surface resistivity function has not been addressed. 

The most critical part of the synthesis procedure is the 
inversion of the (MxN) by (MxN) matrix B. It is desirable to 
use the smallest number (MxN) for computational and accuracy 
reasons. To find the minimum (MxN), a technique is borrowed 
from signal analysis. The scattered field is computed using a 
sufficiently large number of patches. The FFT is then 
performed on this array of field values. The result can be 
translated into the minimum number of patches required to 
accurately represent the field. This is based on the fact that 
plotting a field pattern (function) is basically sampling its 
distribution adequately. The interval between the sampled 
field points at which the field is computed is the sampling 
interval and the reciprocal of this is the sampling frequency. 
Dividing the u and v region into N, and N, segments 
respectively, the corresponding sampling intervals are 2/N, 
and 2/N,. 

Nygquist’s theorem [Ref. 6] requires that at least twice 
the maximum ‘frequency component’ should be used as the 
smallest sampling rate. Utilizing the fact that the FFT is 
symmetric with respect to the half of the sampling frequency 
used, i.e., symmetric to the N,/4 and N,/4 point, a plot of the 
two-dimensional FFT vs the number of segments N, and N, can be 
obtained. For simplicity an equal number of points have. been 
used (N,=N,) and the FFT is performed on the rows and the 


columns of the matrix of the field values to identify the 
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Spatial frequencies in u and v. Figure 24 shows the FFT 
performed on the columns for the field depicted in Figure 12 
using methods (I) and (II). As expected, the frequency content 
is the same in both cases. 

In Figure 25 the spatial frequency content of the field 
in Figure 6 is shown. One can conclude that N, and N, greater 
than 4 1s sufficient. In Figure 26 the field in the ¢=0 plane 
is shown. The dotted line represents the field using only 16 
points (M=N=4) and the solid line the field when M=N=64. 

However for the 5Ax5A sheet with the field shown in 
Figure 12 and the frequency content shown in Figure 27, 6 
segments per direction should be sufficient. In Figures 28 and 
29, the corresponding field plots in the ¢=0 plane show the 
convergence. 

Figure 30 shows the monostatic RCS of a AxA resistive 
Sheet. The spatial frequency content of the field is shown in 
Figure 31, and it reveals that 4 segments are sufficient. The 
corresponding field plots in Figure 32 and 33 illustrate the 
convergence. The same process is applied to the 3Ax3A plate. 
The frequency content is shown in Figure 34 and the 


corresponding convergence is shown in Figures 35 and 36. 
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Figure 25: FFT of the scattered field shown in Figure 6. 
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Figure 26: Bistatic scattered RCS in the #¢=0 plane 
corresponds to Figure 6. The dotted line represent the RCS 
using the synthesized resistivity and M=N=4. 
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Figure 27: FFT in u and v directions of the bistatic 
scattered field in Figure 12. 
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Figure 28: A cut on the u axis corresponding to Figure 12. 
The field (dotted) due to the reconstructed resistivity is 
within a couple of dB of the original. 
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Figure 29: A cut on the u axis corresponding to Figure 12. 
Comparison of original and synthesized patterns for 
M=N=12. 
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Figure 30: Monostatic RCS from a one wavelength square 
resistive sheet. Rs is given by Equation (102). 
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Figure 31: FFT of the monostatic scattered field in Figure 
30:; 
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Figure 32: The monostatic RCS due the reconstructed 
resistivity with 9 patches compared to the original. 
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Figure 33: Using M=N=4, the monostatic RCS due to the 
reconstructed resistivity profile is converged. 
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Figure 34: FFT of the monostatic scattered field from a 
resistive sheet with Lx=Ly=3*lambda. 
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Figure 35: The monostatic scattered field due to the 
reconstructed resistivity (dotted) compared to. the 
original for M=N=4. 
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Figure 36: Using M=N=12 the monostatic scattered field due 
to the reconstructed resistivity is converged. For this 
case Lx=Ly=3*lambda. 
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E. SUMMARY 


In this chapter the synthesis procedure was computer 
Simulated and the corresponding data were presented. Two 
methods were investigated for choosing observation points for 
generating synthesized distributions. In method (I) points 
inside of the visible region were used. This gives good 
results, but it was shown to be very restrictive because the 
number of observation points in the visible region must match 
the number of subdomains on the plate. By allowing complex 
angles it is possible to use observation points in the entire 
region bounded by -isu,vs1l (method (II)). This is much more 
flexible and still gives good results. 

Using Fourier analysis and the sampling theorem it was 
possible to identify the minimum number of patches that can be 
used in the synthesis procedure. In general the larger the 
sheet the more points are necessary, and consequently the 
computing time increases significantly. 

To illustrate the procedure two resistive sheets were 
considered, 1Ax1A and 5Ax5A for the bistatic case and 1Ax1iX 
and 3Ax3A for the monostatic case. The method was tested first 
from the computational point of view. For the resistivity 
distributions considered (constant and linear taper) it was 
found that patch dimensions of approximately 0.5X and 0.25X 
were sufficient in the bistatic and monostatic cases, 
respectively, to obtain good agreement with the original 


scattered field. 
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V. CONCLUSIONS 


A method of RCS synthesis for planar resistive sheets has 
been presented. The synthesis equations were developed 
assuming that the resistive sheet boundary conditions hold at 
all points on the surface. Therefore the equations implicitly 
neglect the effect of the edges on the current. 

The equations were verified by first calculating the RCS 
of the resistive sheet for a known resistivity function and 
then using these values for back solution to obtain a 
synthesized surface resistivity. In principle the two 
resistivity functions should be identical, but because of 
numerical roundoff the agreement is not exact. However, the 
results for a wide range of synthesis parameters were in very 
close agreement with the original resistivity distributions. 

It was also shown that a Fourier decomposition of the 
scattering pattern yields guidelines for choosing the 
subdomain size, or equivalently, the number of far field 
observation points. It was found that a subdomain size of 
about 0.5A for bistatic RCS and 0.25 for monostatic RCS give 
converged results. By convergence is meant that RCS for the 
Synthesized distribution is essentially the same as the RCS 
for the original scattered distribution. 

It 1S important to note that the synthesis is based ona 
particular assumption with regard to the form of the current. 
The resistive sheet boundary condition is equivalent to the 
physical optics approximation in the sense that the current in 
the vicinity of the edge is not disturbed by the presence of 
the edge. The current is only a function of the incident field 
and the surface resistivity (see Equation (4)). Thus using the 
synthesized resistivity distribution in a method of moments 
code could yield a RCS substantially different from that based 
on the assumed current if surface waves exist. Examples are 


shown in Figures 37 and 38. In these cases the surface waves 
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are not as strong as they would be for a surface with more 
highly conducting edges, and therefore the agreement between 
the two curves is relatively close. This problem is overcome 
if the current is treated as an unknown along with the 
resistivity. This approach allows a rigorous solution of the 
problem and is recommended that future work is directed along 


these lines. 


72 





Solid line - MM 


Pe ee ee ee a 


BA tks Se See ee Pets So ND Doo) oe oe EC ON ge A ley eee” 
oe eevee epemuevpnvves anne voevoeesvevevoerense evra snveaneeevvsgePeeevreeeseevenvneQRreeanevoeunvsvnaceseanveererse eee ree eee zr eee ee eee ems ee em eae 
ry . + 


"ee ee ee oe en 6 oe te we ee ee ee ee ee ee wm 


re ee ee a a ee Se 2 ee 2 ee oe ee ed 


Pr a a a a a a a a a a a ar Or er ae © CS Oe ee eS > Ce ee to i ee es 


Pe a er ee 2 er © ce 


‘Sigma/fambda‘2 In dB 


| Ck 
) A | 


Pe ee ee ee | 


arr, , er rs | Mw ea or a re 2 Pe ht eC ee ee i es ee is 
+ 








-~100 -80 # -60 -40 ~20 0 20 40 60 80 100 
Angle theta in deg (phi=0) 


Figure 37: Comparison of the monostatic RCS using MM and 
the approximate method. For this case, Lx=Ly=lambda, 


M=N=12, N,=Nj=12. R, = 3779. 
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Figure 38: Comparison of the monostatic RCS using MM and 
the approximate method. Lx=Ly=lambda, M=N=12, N,=N,=12. R, 
is given by Equation (102). 


74 








APPENDIX. COMPUTER CODES 


This appendix contains the codes developed to computer 
Simulate the synthesis equations. Code field.m calculates the 
scattered field or RCS for each case (monostatic, bistatic, 
method (I), method (II)), given the size of the resistive 
Sheet and the surface resistivity. Code synth.m uses the field 
values calculated from field.m, to synthesize the surface 
resistivity. Code synthld.m uses the surface resistivity 
values calculated by synth.m to plot the RCS in the ¢=0 plane. 
Code fanalysl.m performs the FFT analysis for methods (I) and 
(II), and code fanalys2.m performs the FFT analysis in the u 
and v direction of the scattered field. Function coord.m is 
used in the first two codes to find the coordinated of the 
patches given the size of the sheet and the number of patches. 
Finally, program fit.m finds the closest number of points 
closest to an input value that fit in the unit circle, and is 


used in method (I). 
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% Program field.m 


oe 


* (1) Bistatic case: (sigma_theta_theta, or 


sigma phi phi) /lambda~2 


This program calculates for a rectangular sheet: 


6 or the 

6 scattered field for the general case 
of the 

% incident field, given the values 
Eth i and 

6 Ephi_i and assuming Co=1. 


* (ii) Monostatic case: sigma/lambda*2 


6 The plate can be PEC or resistive. Plots of the 
resistance, 
% as well as, 3-D plots of the field are shown. 


oe 


Works in the direction cosine space. 


oe 


Date: 10 Nov 1994 
By: Nick Faros. 


o\e 


clear 


SESECS INPUT PARAMETERS COCCCS 


choicel=input ('Enter 1 for bistatic or 0 for monostatic: 


76 











choice2=input (’Enter 1 for method (II) or 0 for method (I): 
an ie 


al=input (’Enter size in -x (multiples of wavelength): ’); 


bl=input (‘Enter size in -y (multiples of wavelength): '); 


% Assume wavelength l=1 


L=1% 

k=2*pi/l1; % wavenumber 

a=al*l; &® dimension in x (in meters) 
b=b1*1; % dimension in y (in meters) 


Nx=input (’Enter # of patches orm: oc (N) 2 *-) 3 
Ny=input (‘Enter # of patches in y (M): '); 


dx=a/Nx; - dy=b/Ny; 
xx=-a/2:a/(Nx-1):a/2; =-b/2:b/ (Ny-1) :b/2; 


= xx,yy are the dimension of the plate in meters 


R=input (‘Resistance:0=PEC, l=constant, 2=linear taper, 


3=ramp:’); 

if R== 

cl=input (’Enter constant value: '’); 

end 
tic 
2% 2222S LLLL LEE L LOO LEE EE EEE EE EEE EEE EEE SEES ESE SE EEESESESESS 

22222 FIND COORDINATES OF THE PATCHES %%%%% 

[xo, yo] =coord(a,b,Nx,Ny); % Call function coord.m 
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x1l=reshape (xo,1,Nx*Ny) ; 
yl=reshape (yo,1,Nx*Ny) ; 


COC CCED FIND PATTERN OF THE INPUT RESISTANCE BOCCTCCE CES 


<4 Some cases 


1f R== 
r=zeros (Ny, Nx) ; * Zero everywhere 
elseif R== 
r=cl*ones (Ny, Nx) ; * Constant 
elseif R== 
r=abs (xo) +abs (yo) ; * linear taper in both dimensions 
r=150*rT; 
elseif R== 
r=zeros (Ny, Nx) ; 
ror. L=1.:Nx 
r(:,1)=100*l*ones(Ny,1); % ramp 
end 
end 


ri=reshape(r,1,Nx*Ny) ; 


OoOCO0OC CCC OCOOCOCCOCOC OC OOO OO oOo oOo oOo 88 0000060¢66060606006000,0 0°08 Oo0e00C80.0880 
TTTCOCCTC CD SCOT CCC CCC SC CCCCCESCECCCCECEEEEEEEEEEEEEEEEEEEEELEEESESE SLY 


« Next find coordinates of the center of the patches in 
* direction cosine space. If method (I) is used, run fit.m 


* first, to find nx, ny. If method (II) is used nx=Nx and 


6 ny=Ny; 
if choice2== * Method (IT) 
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nx=Nx; 
ny=Ny; 





ee 


* Method 


elseif choice2==0 


. 
z 


Nu 
Nv 


reshape (u,1,nx*ny) 
reshape (v,1,nx*ny) ; 


ny=5 


end 
[u,v] =coord(2,2,nx,ny) ; 
ul= 
vl= 


a 
f 


sintheta=sqrt (ul.*2 +v1.%2) 
find (sintheta>1) 


I 


° 
z 


sqrt (1-sintheta.%*2) 


atan2 (v1,ul1) 


phi= 


e 
i 


costheta 


% Method (TI) 
costheta (TI) 
u1l(I)=[]; vi(I)=(]; 


if choice2==0 
end 


phi(t)=[]; 


Sintheta(I)=[]; 


1£ choicel==1 


Bistatic 


° 
x) 


e 
i 


") 


translate it into radians 


degrees 


1nN 


input (’Enter theta inc. 


thid 


$ 


thir=thid*pi/180; 


ae 


in degrees 


a11C% 


1 


input (’Enter ph 
phiid*pi/180 


phiid 


* 
i 


phiir 
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coefl=1./(2*r1+377*cos (thir) ); * bistatic 
coef2=1./(2*r1*cos (thir) +377) ; * bistatic 


% u=Sin (theta) *cos (phi) 


6 v=Sin (theta) *sin (phi) 


uil=sin(thir) *cos(phiir) ; 


vi=sin(thir) *sin(phiir) ; 


Sincl=sin(k/2*dx* (u(1,:)+ui))./(k/2*dx* (u(1,:)4+ui)); 
Si=find(isnan(sincl) ); 


Sinc1(S1)=ones(size(S1)); 


Sinc2=sin(k/2*dy* (v(:,1)+vi))./(k/2*dy* (v(:,1)4+wi)); 
S2=find(isnan(sinc2)); 


Sinc2 (S2) =ones(size(S2)); 


Sinc=sinc2*sincl; 
Sinc=reshape (sinc,1,nx*ny) ; 
if choice2== 
sinc(I)=[]; 


end 
choice3=input (’Enter 1 for sigma_theta theta, 2 for 
Sigma_phi_phi, or 3 for the general scattered field: '); 
if choice3==1 * Compute sigma_theta theta 
for p=1:Ny*Nx 


expn=exp (J *k* (x1* (ul (p) tui) +y1* (v1 (p)+4vi))); 
S1(p)=(expn*conj (coefi1’)); 
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end 


%E=dx*dy*cos (thir) *costheta.*sinc.*cos (phi-phiir) 


E=dadx*dy*cos (thir) *costheta.*sinc.*cos (phiir) .*s1; 


RCS1=((377*2*pi*E) .*2) ./pi; 


RCS=10*10g10 (abs (RCS1) ) ; 
RCS=reshape (RCS,Ny,Nx) ; 


elseif choice3==2 * Compute sigma_phi_ phi 


for p=1:Ny*Nx 
expn=exp (j*k* (x1* (ul (p) +ui) +y1* (vl (p) +vi))); 


$2 (p)=(expn*conj (coef2’)); 


end 


E=dax*dy*cos (thir) *sinc.*cos (phi-phiir) .*s2; 
RCS1=(377*2*pi*E) .“2./pi; 


RCS=10*1l0g10 (abs (RCS1) ) ; 
RCS=reshape (RCS ,Ny, Nx) ; 


elseif choice3==3 % Compute field 


Eth i -=input(’Enter magnitude Eth_i: '); 
Ephi_ i =input(’Enter magnitude Ephi_i: ’); 
for p=1:Ny*Nx 
expn=exp (j*k* (x1* (ul (p) tui) +y1* (v1 (p)+vi))); 
s1(p)=(expn*conj (coefl1’ ) ); 


end 
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-*Sl1; 


For p=1:Ny*Nx 
expn=exp (j*k* (x1* (ul (p) tui) +y1* (v1 (p)+vi))); 
S2 (p)=(expn*cony (coef2’)); 

end 


E th=dx*dy*cos (thir) *costheta.*sinc.* (Eth_i*cos(phi-phiir) .* 


S1l+Ephi_ 1*sin(phi-phiir) .*s2) ; 
E phi=-dx*dy*cos (thir) *sinc.* (Eth_i*sin(phi-phiir) .*s1-Ephi _ 
1*cos (phi-phiir) .*s2) ; 


E=sqrt ((abs(E_ th)).*%2+(abs(E_phi)).%*2); 
El=reshape (E,Ny,Nx); 


end 
if choice3== 
save fnb a b E th E_ phi Eth_i Ephi_i r ui vi thir 
phiir ul vl x1 yl xx yy choice3 choice2 % Use in synth.m 
else 
save fnb ab Er ui vi thir phiir ul vl xl yl xx yy 


choice3 choice2 % Use in synth.m 


end 


EEESEEEESES MONOSTATIC CASE $%3%S$%3%3E%S% 
elseif choicel== 


Sincl=sin(k*dx*u(1,:))./(k*dx*u(1,:)); 
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Si=find(isnan(sincl)); 


sinci(S1)=ones (size(S1)); 


sinc2=sin(k*dy*v(:,1))./(k*dy*v(:,1)); 
S2=find(isnan(sinc2) ); 


sinc2 (S2) =ones (size(S2)); 


sinc=sinc2*sincl; 


sinc=reshape (sinc,1,nx*ny) ; 
coef=1./(2*r1+377) ; %¢ monostatic 


for p=1:Ny*Nx 
expn=exp (2*j*k* (x1*ul (p) +y1*v1(p))); 
s(p)=(expn*cony (coef’)); 


end 


E=dx*dy*sinc.*s; 

RCS1= (377*2*pi*E) .*2/pi; 
RCS=10*10g10 (abs (RCS1) ) ; 
RCS=reshape (RCS,Ny,Nx) ; 


if choice2==1 %* Method II 


fe] 


save fnm a b E ul vil x1 yl r xx yy choice2 % Use it in 


synth.m 
elseif choice2==0 %* Method (I) 


E(I)=(]; 
save fom a b E ul vil x1 yl r xx yy choice2 % Use it in 


synth.m 


end 


end 
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save nf E Nx Ny I % Use in fanalys.m 


time=toc 


if choicel== 
if choice2==1 % Method (ITI) 

if choice3==1 | choice3==2 

figure (1) ;mesh(xx, yy, real (r)) ,xlabel (’Size in x’) 


ylabel(’Size in y’),zlabel(’Ohms’) 


figure (2) 

mesh (u,v,RCS) 

xlabel (’u=sin(theta) *cos (phi) ’) 
ylabel ('v=sin(theta) *sin (phi) ’) 
zZlabel (’Sigma/lambda*~2 in dBSm’ ) 


Sine ,.s) 

th=-asin(sint); th=th*180/pi; 

figure (3) 

plot (th, RCS (Ny/2,:)),grid 

xlabel(’Angle theta in deg. at phi=0 plane’) 
ylabel (’Sigma/lambda*2 in dBSm’ ) 


RCS1=reshape (RCS1,Nx,Ny) ; 
mesh(u,v,abs (RCS1) ) 
end 


end 
elseif choicel== 


figure (1) ;mesh(xx,yy,real(r)),xlabel (’Size in 3") 
ylabel('Size in y’),zlabel(’Ohms’ ) 
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figure (2) 

mesh (u,v, RCS) 

xlabel (’u=sin(theta) *cos (phi) ’) 
ylabel (’v=sin(theta) *sin(phi)’) 

- glabel (’Sigma/lambda*2 in dBSm’ ) 


sint=u(1,:)3 

th-asin(sint); th=th*180/pi; 

figure (3) 

plot (th, RCS (Ny/2,:)),grid 

xlabel(’Angle theta in deg. at phi=0 plane’) 
ylabel (’Sigma/lambda*2 in dBSm’) 


RCS1=reshape (RCS1,Nx,Ny) ; 
mesh (u,v,abs (RCS1) ) 

end 

end 

end 

save mf RCS th 
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6 Code synth.m 


« This code reads from the files fn*.mat (output of 

* the code field.m), RCS or field values at specific 

* points, and synthesizes the surface resistivity. 

* This code divides the sheet and DCS into ’big’ 

« patches compared to the field.m code. The number of 
* patches should be choosen such that the center of 

* the big patch, coincides with the center of a small. 
* So, the big patch must contain an odd number of 

* smaller patches. 

* Note that the codes field.m and Synth.m can also 


« have the same number of patches. 


* Works for monostatic and bistatic case. 


clear 


choicel=input (’Enter 1 for bistatic or 0 for monostatic: '); 


if choicel== 
load fnb 
elseif choicel==0 
load fnm 


end 
k=2*pi; 
Nx=input (‘Enter # of patches in x (N): '); 


Ny=input (’Enter # of patches in y (M): '); 


« Nx, Ny are the ’big’ patches. 


tic 
Dx=a/Nx; Dy=b/Ny; 
XX=-a/2:a/ (Nx-1):a/2; YY=-b/2:b/ (Ny-1) :b/2; 
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22222 FIND COORDINATES OF THE PATCHES OF THE SHEET 


[xo, yo]=coord(a,b,Nx,Ny); % Call function coord.m 


x11=reshape (xo,1,Nx*Ny) ; 
y1l=reshape (yo,1,Nx*Ny) ; 


SECCCCS FIND COORDINATES OF THE PATCHES IN DCS 


if choice2== 
nx=Nx; 


ny=Ny; 


[U, V] =coord(2,2,nx,ny); 
U1l=reshape (U,1,nx*ny) ; 
Vi=reshape(V,1,nx*ny) ; 


elseif choice2== 
nx=4 % Nu 


ny=5 6 Nv 


[U,V] =coord(2,2,nx,ny); 
Ul=reshape (U,1,nx*ny) ; 
Vi=reshape (V,1,nx*ny) ; 
sintheta=sqrt (U1.*2 +V1.%2); 


II=find(sintheta>1) ; 
U1 (II)=[]; Vi(II)=[); 


817 


end 


This part finds the index of the selected points 


1:length (U1) 


For 1 


=U le U tL) 


£11 
£21 


=O): 


find (abs (£11) -min(abs(f11) ) 


=V1-V1(1); 


£12 
f22 


find (abs (£12) - 


i 


£22-f£21(11) ; 


dea 


ff= 


find(ff 
ff1==[] 
l14+1 


ffl= 
while 


° 
f 


1i= 


£22-f£21:(11); 


Lf} 


find (ff 


cri= 


end 


af 2: (11. 


index (1) 


end 


o\° 
o\? 
o\? 
ove 
oe 
o\? 
o\? 


ui (index) ; 


UU1 


vi (index) ; 


VVi1 


Ul’ VV1' V1']; 


[UU1’ 
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X1=x1 (index) ; Yl=y1 (index) ; 
ee’ Se YE ya 


sintheta=sqrt (U1.*%2 +V1.%2); 
pphi=atan2 (V1,U1); 


ccostheta=sqrt (1-sintheta.*2); 


if choicel== %* Bistatic 


sinci=sin(k/2*Dx* (U(1,:)+ui)) ./(k/2*Dx* (U(1,:)+ul)); 


Sl=find(isnan(sincl)); 


sinci(S1)=ones(size(S1)); 


sinc2=sin(k/2*Dy* (V(:,1)+vi)) ./(k/2*Dy* (V(:,1)+4+vi)); 


S2=find(isnan(sinc2) ); 


sinc2 (S2) =ones (size(S2) ); 


ssinc=sinc2*sincl; 
ssinc=reshape (ssinc,1,nx*ny) ; 
if choice2== 
ssinc(II)=([]; 
end 
C=([]; 
for l=1:Nx*Ny 
e=exp (j*k* (X1* (U1 (1) tui) +¥1* (V1 (1) +vi))); 
C=[C;e] ; 


end 
if choice3== % Sigma _th th 


EE=E (index) ; 
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E1=EE./ (Dx*Dy*cos (thir) *ccostheta.*ssinc.*cos(pphi-phiinr) ); 
E1=EE./ (Dx*Dy*cos (thir) *ccostheta.*ssinc.*cos(phiinr) ) ; 


Z=find(isnan(E1) ); 
E1(Z)=zeros (size(Z)); 
El=conj (E1’); 
SO=C\EL; 
G=pinv(C)*E1; 
imp=.5*(1./g-377*cos (thir) ); 
impl=reshape (imp, Ny, Nx) ; 


dd=real (imp) ; 


elseif choice3==2 %sigma_phi_ phi 


EE=E (index) ; 

E1=EE./ (Dx*Dy*cos (thir) *ssinc.*cos(pphi-phiir) ); 
Z=find(isnan(E1)); 

Fi (Z)=zeros(size(Z)); 

El=con] (E1’); 

g=C\E1; 

6G=pinv(C) *E1; 

imp=(1/(2*cos (thir) ))*(1./g-377); 
imp1l=reshape (imp, Ny,Nx) ; 


Gd=real (imp) ; 


elseif choice3==3 
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El th=E th./(Dx*Dy*cos (thir) *ccostheta.*ssinc) ; 
Z1=find(isnan(E1_th)); 

Fl th(Z1)=zeros(size(Z1) ); 

El th=conj (E1 th’); 


El phi=E phi./(-Dx*Dy*cos (thir) .*ssinc) ; 
Z2=find(isnan(E1 phi) ); 

El phi (Z2)=zeros (size(Z2)); 

El phi=conj(E1_ phi’); 

El=[El_th;E1_ phi]; 


all=conj (Eth_i*cos(pphi-phiir))’; 


a21=con}j 


( 

al2=conj (Ephi_i*sin(pphi-phiir))’; 
(Eth_i*sin(pphi-phiir) )’; 
( 


a22=conj (-Ephi_i*cos(pphi-phiir) )’; 


Y1l=((al2.*El phi) - (a22.*E1_th))./((a12.*a21-a1l.*a22) ); 
$Y2=(El th-al1l.*Y1)./ai2; 


gli=pinv(C) *Y1i; 
$g2=pinv (C) *Y2; 


imp=.5* (1./g1-377*cos (thir) ) ; 
simp=1/(2*cos (thir) )*(1./g2-377) ; 


impl=reshape (imp, Ny, Nx) ; 


dd=real (imp) ; 


end 
save faros dd Dx Dy k x11 choicel choice3 thir phiir ui 


elseif choicel== % Monostatic 


Sincl=sin(k*Dx*U(1,:))./(kK*Dx*U(1,:)); 
Sl=find(isnan(sinci) ); 


Sinci (S1)=ones(size(S1)); 


sinc2=sin(k*Dy*V(:,1))./(k*Dy*V(:,1)); 
S2=find(isnan(sinc2)); 


Sinc2 (S2) =ones (size(S2)); 


ssinc=sinc2*sincli; 


ssinc=reshape(ssinc,1,nx*ny) ; 


C=[]; 

for l=1:Nx*Ny 
e=exp (2*j*k* (X1*U1 (1) +Y1*V1(1))); 
Ca(Crel? 


end 


EE=E (index) ; 

E1=EE./ (Dx*Dy*ssinc) ; 
Z=find(isnan(E1)) ; 
B1(Z)=zeros(size(Z)); 
El=conj (E1’); 

g=C\E1l; 

6g=pinv(C) *E1; 
imp=.5*(1./9-377) ; 
imp1l=reshape (imp, Ny, Nx) ; 


dd=real (imp); 


° 


save faros dd Dx Dy k x11 choicel % Use in synthid.m 
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time=toc 


end 


figure 

mesh (XX, YY, real (imp1) ) 

hold on 

mesh (xx, yy,TY) 

hold off 

xlabel(’x (wavelengths) ’) 
ylabel(’y (wavelengths) ’”) 
zlabel(’Rs and Ros in Ohms’ ) 


rerror=max (max (abs (imp1) -r) ) /max(max(r) ) 
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% Code synthid.m 
* This code reads the surface resistivity from the file 
faros.dat 
* created by field.m and plots the scattered RCS in the 
principal 


6 plane phi=0 (u axis). 


clear 


load faros 


Sinth=0: (1/49) :1; % Sampling the sin (or the dcs) 
costh=sqrt (1-sinth.%2) ; 


if choicel==1 % Bistatic 


1£ choice3== * Sigma_th_th 


G=1./ (2*dd+377*cos (thir) ); 
fOr p=1:50 
expn=exp (j}*k* (x11* (sinth(p)+ui))); 
S(p)=(expn*conj (d)); 
end 
sincl=sin(k/2*Dx* (sinth+ui) )./(k/2*Dx* (sinth+ui) ); 
S1=find(isnan(sincl) ); 


Sinc1(S1)=ones(size(S1)); 


RCS= (377*2*pi*Dx*Dy*cos (thir) *costh*cos (phiir) .*sincl.*s) .*2 


/pi; 


elseif choice3==2 % Sigma phi phi 
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d=1./(2*dd*cos (thir) +377) ; 

for p=1:50 
expn=exp (j*k* (x11* (sinth(p)+ul))); 
s(p)=(expn*cony (d)); 


end 


sinci=sin(k/2*Dx* (sinth+ui)) ./(k/2*Dx* (sinth+ui) ); 
Si=find(isnan(sinci)); 
sinci (S1) =ones (size(S1)); 


RCS= (377*2*pi*Dx*Dy*cos (thir) *cos (phiir) *sincl.*s) .*2/pi; 


elseif choice3== 


brake 


end 


elseif choicei==0  Monostatic 


d=1./(2*dd+377); 
For p=1:50 
expn=exp (2*j}*k* (x11*sinth(p))); 
s(p)=(expn*cony (d)); 
end 


sincl=sin(k*Dx*sinth) ./(k*Dx*sinth) ; 


S1l=find(isnan(sincl)); 


sincl1(S1)=ones(size(S1)); 
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RCS= (377*2*pi*Dx*Dy*sincl.*s) .*2/pi; 


end 


RCS1=fliplr (RCS) ; 
RCS=[RCS1 RCS]; 


theta=asin(sinth) ; 
thetal=fliplr (theta) ; 
theta2=[-thetal theta]; 
theta3=theta2*180/pi; 
Sinthetal=fliplr(sinth) ; 
Sintheta2=[-sinthetal sinth]; 


figure (1) 
plot (theta3,10*1lo0g10(abs(RCS)),’r’) 
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% Code fanalysl.m 
This code calculates the FFT of the field 
> of methods (I) and (II) 





ol? 


Ne 


clear 

load nf 

% Method (TI) 
x=length(E) ; 

Yl=fft (E) ; 

MM= [Nx Ny] ;MMM=min (M™) ; 
£1=MMM/2* (0:x/2-1) /x; 


% Method (ITI) 

E2=E; 

E2(I)=(1; 
x2=length(EZ2) ; 

Y=fft (E2) ; 

f2=MMM/2* (0:x2/2-1) /x2; 


figure, subplot (211) 

plot (£2,abs (Y(1:x2/2) )) 

saxis([0 3 0 1500]); 

xlabel (’Nu or Nv’);ylabel(’abs (FFT) ’) 


subplot (212) 
plot (f1, (abs (Y1(1:x/2)))); %taxis([0 3 0 1500]); 


xlabel ('Nu or Nv’) ;ylabel (’abs (FFT) ’) 
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6 Code fanalys2., 
6 This code calculates the FFT of the field 


* in u and v directions 


clear 
load nf 


* Method (I) columns 
SE=(377*2*pit*E) .*2/pi; 
x=length (E) ; 
Yi=abs (fft (E)); 

MM= {Nx Ny] ;MMM=min (MM) ; 
f1=MMM/2* (0:x/2-1)/x; 


* Method (II) -rows 
k2=reshape (E,Ny,Ny) ; 
E2=conj (E’); 

E2=reshape (E2,1,Nx*Ny) ; 
Y2=abs (fft (E2)); 


figure, subplot (211) 

plot (fi, Yit12ex/2) ) 

axis({[0 5 0 .8)); 

xlabel(’ Nv’);ylabel(’abs(FFT) ’) 


subplot (212) 
PLOU(£1, Yi x72) )% 


axis([0 5 0 .8)); 


xlabel (’Nu’) ;ylabel (’abs (FFT) ’) 
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function [xo, yo] =coord(a,b,Nx,Ny) 

% Given the dimensions (a) in x, (b) in y of a 

% rectangle and the number of segments (Nx) in x 
% and (Ny) in y, the coordinates of center of 

* each patch is returned. 

% Point 0 of the coordinate system is at the 


6 center of the rectangle. 


* Date: 1 Nov 1994 
* By: Nick Faros 


dx=a/Nx; 
dy=b/Ny; 


if rem(Nx,2)==0 ¢ 1.e. Nx is an even number 


for n=1:Nx 


for 1=1:Ny 
xo (i,n) =- (Nx/2- (n-1) -0.5) *dx; 
end 
end 
else % if Nx is an odd number 


for n=1:Nx 


fOr i=1:Ny 
xo (i,n)=- ((Nx-1) /2- (n-1) -0.5) *dx-dx/2; 
end 
end 


end 

if rem(Ny,2)==0 6 i.e. Ny is an even number 
for n=1:Nx 
for i=1:Ny 


yo(i,n)= (Ny/2-(i-1)-0.5) *ay; 


end 
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end 


else %* 1f Nx is an odd number 


Lor =i 2Nx 
Lor 2=1:Ny 
yo(i,n)= ((Ny+1)/2- (i-1)-0.5)*dy-dy/2; 
end 

end 


end 


100 








Code fit.m 
6 This program fits M points in the unit circle. 


oe 


oe 


The spacing is uniform in one or both directions. 


6 Date: 1 Nov 1994 

S By: Nick Faros. 

clear 

M=input (‘Enter number of points in the unit circle: '); 


NnX=2; ny=2; %¢ Start at those values. 


Aax=2/nx; 


dy=2/ny; 


% Find the coordinates of the center 


% of the patches. 


[ul,vi] =coord(2,2,nx,ny)j; 
ul=reshape (ul,1,nx*ny) ; 
vi=reshape (v1,1,nx*ny) ; 
sintheta=sgqrt (ul.*2 +v1.%2); 
T=find(sintheta>s1) ; 
vi(I)=[];ul(1)=[]; 
L=length (v1) ; 

k=1; 


while L<M 
if rem(k,2)== 
nx=nx+1; 
[ul ,v1]=coord(2,2,nx,ny); 
ul=reshape (u1l,1,nx*ny) ; 
vi=reshape (v1,1,nx*ny) ; 
sintheta=sqrt(u1l.*2 +v1.%2); 
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T=find(sintheta>1) ; 
v1 (I) =f[} rab es i tl 
L=length (v1) ; 


elseif rem(k,2)~=0 
ny=ny+1; 

[ul,v1] =coord(2,2,nx,ny); 
ul=reshape (ul,1,nx*ny) ; 
vi=reshape(v1,1,nx*ny) ; 
Sintheta=sqrt (ul.*2 +v1.%2); 
I=find(sinthetasl1) ; 
Vitlyeatl ult tyat)> 
L=length (v1); 

end 

k=k+1; 

end 

[nx ny L] 

plot (ul,v1,’*'),xlabel(’u’),ylabel(’v’) 


title(’Points inside the unit circle’) 
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